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Abstract 

In this article, we assume that the nonet scalar mesons below 1 GeV are the two-quark- 
tetraquark mixed states and study their masses and pole residues using the QCD sum rules. 

In calculation, we take into account the vacuum condensates up to dimension 10 and the 
0{as) corrections to the perturbative terms in the operator product expansion. We determine 
the mixing angles, which indicate the two-quark components are much larger than 50%, then 
obtain the masses and pole residues of the nonet scalar mesons. 
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1 Introduction 

There are many scalar mesons below 2 GeV, which cannot be accommodated in one qq nonet, some 
are supposed to be glueballs, molecular states and tetraquark states [mil IIS]. In the scenario of 
molecular states, the scalar states below 1 GeV are taken as loosely bound mesonic molecular states 
d, or dynamical generated resonances [7]. On the other hand, in the scenario of tetraquark states, 
if we suppose the dynamics dominates the scalar mesons below and above 1 GeV are different, 
there maybe exist two scalar nonets below 1.7 GeV n n u. The strong attractions between 
the scalar diquarks and anti-diquarks in relative S'-wave maybe result in a nonet tetraquark states 
manifest below 1 GeV, while the conventional quark-antiquark nonet mesons have masses about 
(1.2 —1.6) GeV. The well established ^Pi and ^P 2 quark-antiquark nonets lie in the same region. In 
2013, Weinberg explored the tetraquark states in the large-Vc limit and observed that the existence 
of light tetraquark states is consistent with large-Vc QCD |8] . We usually take the lowest scalar 
nonet mesons {/o/cr(500), ao(980),/to(800),/o(980)} to be the tetraquark states, and assign the 
higher scalar nonet mesons {/o(1370), ao(1450), Kq (1430),/o(1500)} to be the conventional ^Pq 
quark-antiquark states miiiin]- 

There maybe exist some mixing between the two scalar nonet mesons, for example, in the chiral 
theory [10]. In the naive quark model, for /o(980) = ss, the strong decay /o(980) —^ tttt is Okubo- 
Zweig-Iizuka forbidden; for aQ(980) = , the radiative decay ^(1020) —>■ aQ(980)7 is both 

Okubo-Zweig-Iizuka forbidden and isospin violated. From the Review of Particle Physics, we can 
see that the process /o(980) —>■ tttt dominates the decays of the /o(980) and the branching fractions 
Br (((>(1020) ^ oj] (980)7) = (7.6 ± 0.6) x 10-^ Br ( ((.(1020) ^ /o(980)7) = (3.22 ±0.19) x lO’^ 
[ij. The naive quark model cannot account for the experimental data even qualitatively, we have 
to introduce some tetraquark constituents, such as and jf do not want to 

turn on the instanton effects m- 

We can use QCD sum rules to study the two-quark and tetraquark states. QCD sum rules 
provides a powerful theoretical tool in studying the hadronic properties, and has been applied 
extensively to study the masses, decay constants, hadronic form-factors, coupling constants, etc 
[mils]. There have been several works on the light tetraquark states using the QCD sum rules 
[HdSKisKnKiHKiiiiiniiiiiiii]- inRefs.[iiiin], the scalar nonet mesons below 1 GeV are taken 
to be the tetraquark states consist of scalar diquark pairs and studied with the QCD sum rules 
by carrying out the operator product expansion up to the vacuum condensates of dimension 6. In 
Ref. [18] , Lee carries out the operator product expansion by including the vacuum condensates up 
to dimension 8, and observes no evidence of the couplings of the tetraquark currents to the light 
scalar nonet mesons. In Ref. m, Chen, Hosaka and Zhu study the light scalar tetraquark states 
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with the QCD sum rules in a systematic way. In Ref.[5D], Sugiyama et al study the non-singlet 
scalar mesons ao(980) and ko(800) as the two-quark-tetraquark mixed states with the QCD sum 
rules, and observe that the tetraquark currents predict lower masses than the two-quark currents, 
and the tetraquark states occupy about (70 — 90)% of the lowest mass states. 

In this article, we assume that the scalar nonet mesons below 1 GeV are the two-quark- 
tetraquark mixed states and study their properties with the QCD sum rules in a systematic way 
by taking into account the vacuum condensates up to dimension 10 and the 0(as) corrections to 
the dimension zero terms in the QCD spectral densities in the operator product expansion. 

The article is arranged as follows: we derive the QCD sum rules for the scalar nonet mesons in 
Sect.2; in Sect.3, we present the numerical results and discussions; and Sect.4 is reserved for our 
conclusions. 


2 The scalar nonet mesons with the QCD Sum Rules 


In the scenario of conventional two-quark states, the structures of the scalar nonet mesons in the 
ideal mixing limit can be symbolically written as 


/o(500) = 


dd 


’ 

Oq (980) = du, 
Kq (800) = us, 


/o(980) = ss, 

n UU — dd 

a°(980) = 


V2 ’ 

«:Q(800) = ds, «:q( 800) = sd 


Oq (980) = ud, 


Kq (800) = su. 


( 1 ) 


In the scenario of tetraquark states, the structures of the scalar nonet mesons in the ideal 
mixing limit can be symbolically written as 


fo{500) = udud , /o(980) = 


Oq (980) = dsus, 
(800) = udds , 


a[J(980) = 


usus -\- dsds 
usus — dsds 


Kq(800) = udus , «:q(800) = usud 


a(j" (980) = usds , 


Kq (800) = dsud. 


( 2 ) 


If we take the diquarks and antidiquarks as the basic constituents, the two isoscalar states udud 
and mix ideally, the ss degenerates with the isovector states ssdu, ss and 

ssud naturally. The mass spectrum is inverted compare to the traditional qq mesons. The lightest 
state is the non-strange isosinglet, the heaviest states are the degenerate isosinglet and isovector 
states with hidden ss pairs, the four strange states lie in between. 

In this article, we take the scalar nonet mesons to be the two-quark-tetraquark mixed states, 
and write down the two-point correlation functions ns(p), 


ns(p") = iJd^xe^P--{0\T[js{x)JsHQ)}\0), (3) 

Js{x) = cos9sJg{x)+smesJs{x), (4) 
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where S = /o(980), aQ(980), «;J(800), /o(500), and 


(4 

'-'/o(980) 

(a;) 

•^/oigso) 

(x) 

00 

00 

0 

(x) 

(2 

-'a°(980) 

(x) 

74 

K^SOO) 

(x) 

/2 

(x) 

t4 

■^fo{500) 

(x) 

"'/o(500) 

(x) 


^ijk ^imn 


V2 

M., 

3^2 


^ijk ^imn 


{uj {x)C-f5Skix) Umix)j5Cs'^{x) + dj {x)Cj5Skix) dmix)-f5C s"^ (x)} , 
s(x)s(a;), (5) 


^ {uJ{x)C-f5Skix)Umix)j5Cs'^{x) - dj {x)Cj5Skix) dm{x)-/5Cs^{x)} 
(ss) u{x)u{x) — d(x)d(x) 


6 72 

{x)C-f 5 dk{x) Sm{x)^bCd^{x), 


{m) -I . , X 

- , 

6 

^zjk^tmn (^x)Cj 5 dkix) Um{x)-f 5 Cd^ (x) , 
(qq) u{x)u{x) + d(x)d(x) 


372 


72 


( 6 ) 

(7) 

( 8 ) 


the currents Js{x) and J^{x) are tetraquark and two-quark operators, respectively, and couple 
potentially to the tetraquark and two-quark components of the scalar nonet mesons, respectively, 
the 9s are the mixing angles. In the currents Jg(x), the i, j, k, ... are color indices and C 
is the charge conjugation matrix, the {x)C^^dk{x), {x)C^ 5 Sk{x), e^^^dj{x)C^ 5 Sk{x) 

represent the scalar diquarks in color anti-triplet, the corresponding antidiquarks can be obtained 
by charge conjugation. The one-gluon exchange force and the instanton induced force can result 
in significant attractions between the quarks in the scalar diquark channels [31 [33] . 

In the following, we perform Fierz re-arrangement to the currents and >77(980) both in 

the color and Dirac-spinor spaces to obtain the result, 


J' 


/o(980) 


If _ UU + dd _ ui^5U + di'f^d 
— ■( —ss -;=-h SlJsS- 


72 


72 


- u uj^u + djf,d _ u'yfj,^5U + d'yfj,^5d 

sq'^s—-— -j= — -sy'^yss — 


72 


72 


1 _ ua^'^u + da^'^d su us + sd ds si^^u uij^s + sij^d di^^s 
^ s - i ^^ uu ' y^s + sjf ^ dd ' j^s ^ 57^75^117 ^755 -I- 57^75^(17^755 


72 

1 sa^vuuu^''s -\- su^^dda^'^s 

2 72 


72 


(9) 




ag(980) 


UU — dd 


uij5U-di^5d _ u^^u-d'y^d _ u'y^'y^u - d'y^.'y^d 
— ss -;=-h 51755 - — -S 7'^5 - - - — - - - 57^^755 - - - - - 


72 


72 


72 


72 


1 _ ua^'^u — da^’^d .su us — sd ds si'y^u uij^s — sij^d di^^s 

+ -j^sa^yS - 7 =-h 


72 


72 


72 


s ^ f ^ nuj^s - s '^^ d d ' y^s ^ 57^75^1(7 ^755 - 57^75(^(17^755 


72 


72 


1 sa^vuuu^’' s — sa^^dda^" s 


72 


( 10 ) 


some components couple potentially to the meson pairs tttt, KK, lyn, the strong decays /o(980) —> 
TTTT, KK and ao(980) —>■ rjir, KK are Okubo-Zweig-Iizuka super-allowed, which can also be used 
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to study the radiative decays (/)(1020) —>■ /o(980)7 and (^(1020) —?> aQ(980)7 through the virtual 
KK loops. So it is reasonable to assume that the nonet scalar mesons below 1 GeV have some 
tetraquark constituents. 

The tetraquark operator Jg{x) contains a hidden qq component with q = u, d or s. If we con¬ 
tract the corresponding quark pair in the currents Jgix) and substitute it by the quark condensate 
B then 


'^/o(980) (^) 


'^/o(980)(^) ’ 


'^aO(980) (^) 

-1 

'^og(980)(^) ’ 


•^k+(800) 

-1 

'^k+(800)(^) ’ 


’^foi500) (^) 

-1 

'^/o(500)(®) • 

( 11 ) 


The contracted parts appear as the normalization factors 

currents (a;), ^«o( 800 )(^) and respectively. 

We insert a complete set of intermediate states with the same quantum numbers as the current 
operators Jsix) satisfying the unitarity principle into the correlation functions n 5 (p^) to obtain 
the hadronic representation mis]. After isolating the ground state contributions from the pole 
terms of the scalar nonet mesons, we get the result, 

Us{p^) = ^^ + --- , (12) 

rrig — 

where we have used the definitions (0| Js(0)|5') = As for the pole residues. 

The correlation functions can be re-written as 


ns(p^) = cos^ 9Ug^{p'^) + sin0cos0n|^(p^) -|- sin 0 cos 0 ) -|- sin^ 0n|^(p^), 

nr(p') = ^Jd^xe^P■^{0\T{J^ix)Jf{0)}\0), (13) 


where m,n = 2,4. We can prove that Wg‘^{p'^) = ng'"(p^) with the replacements x —>■ —x and 
p —> —p for m ^ n. 

In the following, we briefly outline the operator product expansion for the correlation functions 
n™”(p^) in perturbative QCD. Firstly, we contract the u, d and s quark fields in the correlation 


^ For example, 


r4 

-'/o(980) 


~vr~ 

V2 

^ijk^imn 


{- [G75]„^ [ 75 C];,, NG],. 




^jm^Xct /_ V ^nk^Tj3 _ 


- {uu) 


12 


SS - [<^ 75 ] c /3 [75 C];^^ 




12 


{dd) 


{uu) + {dd) 

247/2 


Tr {[C75] [75 C]^ '^ss = -^ss = J 


3\/2 


t2 

■'/o(980) ^ 


^nk^rS _ 1 

- —SS > 

12 / 


where the q, /9, A and r are Dirac spinor indexes. 
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functions 11™" (p^) with Wick theorem, and obtain the results: 

{Tr [j 5 Skk'{xh 5 CUjj,{x)C] Tr [y 5 Sn'n{-x)-f 5 CU^,^{-x)C] 

+Tr [j 5 Skk'ix)j 5 CDfj,(x)C] Tr [y 5 Sn'ni-x)-f 5 CD'^,^{-x)C] } , 

n«o( 8 oo)(p^) = j d^x e^P-^ 

Tr [-f^Dkk'{x)-i8CUjj,{x)C] Tr [-fzDn-n{-x)-i^CS'^,^{-x)C] , 

n/o(500)(p^) = j d^x dP-^ 

Tr [-f^Dkk-{x)-i8CUjj,{x)C] Tr [-fzDri-n{-x)-i^CU'^,^{-x)C] , (14) 

n/o( 980 )(P^) = y d'^a; e*P'^Tr[S'jfc(a;)S'fcj(-a;)] 

+ J d*x dP-^Ty[Ska{x)San{-x)an , 

n^^( 980 )(/) = dP-^ {Tr [U,k{x)Ukj{-x)] + Tr [D,k{x)Dkj{-x)]} 

+ J d^x dP-- 

(Tr [U,aix)Ua„,i-x)aP‘'] + Tr [i^,„(x)7^am(-x)a^"]} , 
nf„(800)(P^) = ld^xdP-^Tr[U,k{x)Sk,{-x)] 

+ ^d^’^d^^{qna^,qk)i J d*x dP-^TT[U,a{x)Sa^{-x)aP''] , 

KdooM) = * / d^x (Tr [U,kix)Ukji-x)] + Tr [D,kix)Dk,i-x)]} 

(Tr [Uja{x)UaU-x)<Jp''] + Tr [i^,„(x)i^am(-x)a^"]} , (15) 

n/o(980)(^' ) = n^g(ggQ)(p ), 

n^o(980)(p") = K(9SoM)^ 

K(SOO)iP^) = nf„(goo)(/), 

n^o(500)(p") = n4o(500)b^), (16) 
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n 


22 

/o(980) 


(P^) 


18 


■i I d'^x e'P''^Tr[Sjk{x)Skj{-x)\ 


nao(980)(P^) 
nKo(800)(P ) 
n^o(500)(P^) 


e^^'"{^AU,k{x)Ukj{-x)]+TT[D,k{x)Dkj{-x)]} , 

J e*P-"Tr[[/,fc(x) 5 fe,(-x)] , 

+ Tr [Djk{x)Dkj{-x)]} , (17) 


where 
[/ij {x) 


Dij {x) 


Sij {x) 


iSjj ^ _ dijiriq _ Sij{qq) 
2 tt'^x‘^ Att'^x'^ 12 

32'!t‘^x‘^ 

Uijiyx) , 


'idij ^ SijTTls diji^Ss') 

27r^x‘* 47r^x^ 12 


327r^x^ 


iSjj ^mq{qq) _ SjjX^{qgsuGq) 

48 ~ 192 

- H-, 

iSij ^ms(ss) 5ijx‘^{sgsaGs) 
48 192 

“ 4“ * * * 5 

O 


iSjjX^ ^mq{qgsaGq) 

1152 


iSijx'^ ^ms{sgs<jGs) 

1152 


(18) 


where q = u, d[T3]. We take the assumption of vacuum saturation for the higher dimension vacuum 
condensates and factorize the higher dimension vacuum condensates into lower dimension vacuum 
condensates [12], for example, {qqqq) ^ {qq){qq), {ggdgsO'Gq) ^ {qq){qgsaGq), where q = u,d,s. 
Factorization works well in large Nc limit, in reality, Nc = 3, some (not much) ambiguities maybe 
originate from the vacuum saturation assumption. 

In Fig.l, we show the Feynman diagrams containing the qq annihilations accounting for the 
mixing of different Fock states. The quark-pair annihilations are substituted by the condensates 
{qq){q'q') and {qq){q'gscGq') as there are normalization factors {qq) in the interpolating currents 
Jj (x). The perturbative part of the quark-pair annihilations must disappear as only the terms (qq) 
and {djU^'^qi) in the full quark propagators Uij{x), Dij{x) and Sij{x) survive in the limit x —>■ 0, 
where q = u,d, s. 

In Eq.(18), we retain the terms {djO^^qi) and {sjU^^Si) come from the Fierz re-arrangement of 
the {qidj) and {siSj) to absorb the gluons emitted from other quark lines to form {dj 9 sG'^^t'^„afj_,yqi) 
and {sjgsG'^pt'^^a^vSi) to extract the mixed condensates (dgsO'Gq) and (sgsaGs). Some terms 
involving the mixed condensates (dgsO'Gq) and (sggoGs) appear and play an important role in 
the QCD sum rules, see the second Feynman diagram shown in Fig.l and the first two Feynman 
diagrams shown in Fig.2. 

Then we compute the integrals in the coordinate space to obtain the correlation functions 
ns(p^), therefore the QCD spectral densities ps{s) at the quark level through the dispersion 
relation, 


Ps{s) 


Imn(s) 

TT 


In this article, we approximate the continuum contributions by 


(19) 


dsps{s) exp(^-^'j , ( 20 ) 

which contain both perturbative and non-perturbative contributions, we use the Sg to denote the 
continuum threshold parameters. For the conventional two-quark scalar mesons, only perturbative 
contributions survive in such integrals, see Eqs.(26-27), Eq.(30) and Eq.(33). 
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Figure 1: The Feynman diagrams contribute to the condensates {qq){q'q') and {qq){q'gauGq') in 
the correlation functions where q,q' = u,d,s and S = /o(980), ao(980), ko(800), /o(500), 

the large • denotes the normalization factors (qq) in the currents J|(0). Other diagrams obtained 
by interchanging of the quark lines are implied. 






Figure 2: The Feynman diagrams contribute to the condensates {qgsoGq) and {qq){q'gsuGq') in 
the correlation functions where q,q' = u,d,s and S = /o(980), ao(980), ko(800), /o(500). 

Other diagrams obtained by interchanging of the quark lines are implied. 
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In this article, we carry out the operator product expansion by including the vacuum con¬ 
densates up to dimension 10. The condensates {g^GGG), { °‘^^^ ){qgsa-Gq) have the 

dimensions 6, 8, 9 respectively, but they are the vacuum expectations of the operators of the or¬ 
der 0 (q!s), 0{o?J‘^) respectively, their values are very small and discarded. We take the 

truncations n < 10 and fc < 1, the operators of the orders 0{a’^) with fc > 1 are discarded. Further¬ 
more, we take into account the 0{as) corrections to the perturbative terms, which were calculated 
recently [52]. As there are normalization factors {qq)'^ in the correlation functions H^{p) , we 
count those perturbative terms as of the order (qq)^, and truncate the operator product expansion 
to the order (qq)'^{q'q'), where q, q' = u, d, s. 

Once the analytical QCD spectral densities are obtained, then we can take the quark-hadron 
duality below the continuum thresholds Sg and perform the Borel transformation with respect to 
the variable = —p^, finally we obtain the QCD sum rules, 



psis) = cos^ 9s Ps^(s)+2 sinOs cos 9s p‘^(s) + sin^ 9s p'^si^) ^ ( 22 ) 


P/o/ao(980) 


Pfo{980) 

Pao{980) 

22 

P/o(980) 

22 

Pao(980) 


as f 57 p? 

1 + — —+ 21og^ 

TT \ 5 s 


{mq - 2ms){qq) + {ms - 2m,)(ss) 


614407r6 

(3ms - mq){qgsaGq) + {3mq - ms){sgsaGs) 


1927r4 

, {q<i){ss) 
s H- , s 


1927r4 127r2 

{qq){sgsaGs) + {ss){qgsaGq) {qgsaGq){sgs(jGs) 


247r2 967r2 

{2mq - ms){qq){ss)'^ + {2ms - mq){ss){qq)'^ 


-6{s) 


6{s) 


ISSGtt"^ tt 


^ asGG ^ 


ms{qq) + mq{ss) asGG mq{qq) + ms{ss) .agGG 


727r2 


-) + 


1927r2 


+ )S{s ), 

216 TT 


144 1 TT^ 


■s + {— - )d{s) + 24ms(ss)(5(s) 


(s-s)2 


288 1 7r2 




-p { ‘^‘'^^ )6{s) + 2Amq{qq)5{s) 


144 1 7r2 


TT 

, asGG 

TT 


)5{s) + 24ms(ss)(5(s) 


(s-s)2 


288 7r2 


^.g -p { ^'‘^^ )5{s) + 2Amq{qq)5{s) 


{qq){qgs(jGq) 

967r2 

{ss){sgsCrGs) 

1927r2 


(23) 

(24) 

(25) 

(26) 
(27) 




























44 

Pkq{800) 


Pkq{800) 


Pk,q{800) 


f T , {ms-2mq){ss) - {■mg + 2ms){qq) ^ 

6l3ito5i' + TU + ^‘°*Tj| +- 383^3-* 

3{ms + mq){qgsaGq) + (3mg - ms){sgsCrGs) (qq)'^ + {qq){ss) 

3847r‘* ^ 2471^ 

‘^{m){qgsCrGq) + {qq){sgsaGs) + {ss){qgsaGq) 


4877^ 

{qgsaGq)'^ + {qgs(jGq){sgsCrGs) 

1927r2 

(2ms - mq){qq)^ + {Aniq - ms){ss){qq)'^ 


Sis) 


18 




UsGG (ms — 2mq){ss) — (niq + 2ms){qq) .asGG 


-) + 


153671^ TT 38477^ 

i2mq + ms){qq) + mq{ss) ^agGG^ 




5767r2 


-(- 


-) 


+tL [{gg)^ + {gg){ss)] , 


432 

M 

288 1 7r2 


(28) 


3 Q GG 

—2 s + {— -)<5(s) + 4(mg + 2ms){qq)S{s) + 4(ms + 2mq)(ss)6{s) 


{gq){ggsCrGq) 

^ 192772 

I-^s + + 4(mg + 2ms){qq)dis) + 4(ms + 2mq)(ss)5(s)l 

ZOO TT^ TT J 


(29) 

(30) 


Pfo{500) 


Pfo{500) 
22 

^/o(500) 


6 I 44 O 776 



57 fj.^ 

— + 2 log — 
0 s 


mqjqq) 2 (gg)^ mqjqgsaGq) 

9077^ 12772 ^ 4877^ 


{qq){qgsaGq) (qgsaGq)'^ 


12772 


96772 


<5(s) - 


2mq{qq)^ 


Sis) 


ttgGG. 5mq{qq) agGG 


153677 ^ 77 


-)- 


a,GG 


288772 


) + :^(g-g)^( 

TT 216 TT 


)<5(s) 


(g-g)^ 


144 1 772 


+ { °^‘'^^ )6is) + 24m,(gg)(5(s) 


{qg)iqgscrGq) 

96772 


(g-g)^ 


144 1 772 


+ { °^‘'^^ )6is) + 24:mq{qq)Sis) 


(31) 

(32) 

(33) 


We differentiate Eq.(21) with respect to — then eliminate the pole residues As, and obtain 
the QCD sum rules for the masses, 


2 /o ® Ss Psis) exp (-^) 

ms = - 


/o^ dspsis) exp(-j^) 


(34) 


3 Numerical results and discussions 

In calculation, the input parameters are taken to be the standard values (ss) = (0.8 ± 0.2){qq), 
{sgsCrGs) = mg(ss), {qgsCrGq) = mliqq), mg = (0.8 ± 0.2) GeV^, {uu) = (dd) = (qq) = —(0.24± 
O.OlGeV)^, {^^^^) = (0.33 GeV)"*, m„ = m^ = 6MeV and m^ = 140MeV at the energy scale 
/i = 1 GeV [13 [131 m] . The values m^ = m^ = 6 MeV can also be obtained from the Gell-Mann- 
Oakes-Renner relation at the energy scale /i = 1 GeV in the isospin limit. 
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ms (MeV) 

Fs (MeV) 

ms + rs/2 (MeV) 

ws-rs/2(MeV) 

/o(980) 

990 ± 20 

40 - 100 

1025'= 


/o(1500) 

1504 ± 6 

109 ±7 


1450" 

ao(980) 

980 ± 20 

50 - 100 

1018'= 


ao(1450) 

1474 ±19 

265 ± 13 


1342" 

ko(800) 

682 ± 29 

547 ± 24 

956" 


7s:^(1430) 

1425 ± 50 

270 ± 80 


1290" 

/o(500) 

400 - 550 

400 - 700 

750" 


/o(1370) 

1200 - 1500 

200 - 500 


1250* 


Table 1: The Breit-Wigner masses and widths of the scalar mesons from the Particle Data Group, 
where the superscript c denotes the central values, and the superscript * denotes that we have taken 
the low bound of the width of the /o(1370). 


Firstly, let us set the mixing angles 9s in the QCD spectral densities ps{s) in Eq.(22) to be 
zero, then the scalar nonet mesons are pure tetraquark states. The perturbative QCD spectral 
densities are proportional to it is difficult to satisfy the pole dominance condition PC > 50% 
if the continuum threshold parameters Sg are not large enough and the Borel parameters are 
not small enough, where the pole contribution (PC) is defined by 


lo" dspsjs) exp(-^) 

ir^^Psis) exp(-^) ■ 


(35) 


For Sg, it is reasonable to take any values satisfying the relation, Wgr + < wist — , 

where the gr and 1st denote the ground state and the first excited state (or the higher resonant 
state) respectively. The \/s 5 lies between the two Breit-Wigner resonances, if we parameterize 
the scalar mesons with the Breit-Wigner masses and widths. More explicitly. 


W/o(980) + 
Wao(980) + 
Wko( 800 ) + 
W/o(500) + 


r/o(980) 

2 

Too (980) 
2 

rKo(800) 

2 

r/o(500) 

2 


< 

< 

< 

< 



0 

/o(980) 


0 

ao (980) 


0 

ko(800) 


0 

/o(500) 


— ’^/o(1500) ~ 

< Wap (1450) — 

< Wi^;'.(1430) - 

< W/p(i37o) - 


r/Q(1500) 

2 

rao(1450) 

2 ’ 
rK*(1430) 
2 

^/o(1370) 

2 


(36) 


In Table 1, we show the Breit-Wigner masses and widths of the scalar mesons from the Particle 
Data Group explicitly [1]. Based on the values in Table 1, we can choose the largest continuum 
threshold parameters = 1.9 GeV", sO^(ggo) = 1.8 GeV^, = 1.7 GeV^ and 3°^^^^^^ = 

1.6 GeV^ tentatively to take into account all the ground state contributions and avoid the possible 
contaminations from the higher resonances /o(1370), ao(1450), iFo(1430) and /o(1500). 

In Fig.3, we plot the masses of the scalar mesons as pure tetraquark states with variations 
of the Borel parameter where the central values of other parameters are taken. From the 
figure, we can see that if we exclude the contributions of the condensates {qq){q'ggCrGq') with 
q, q' = M, d, s, the predicted masses ms increase monotonously and quickly with increase of the 
Borel parameters at the value < 0.9 GeV^, then increase slowly and reach the values 
”^/o( 980 ) = 1.06 GeV, Wao( 980 ) = 1.03 GeV, w,.,p( 8 oo) = 0.99 GeV, w/o( 5 oo) = 0.96 GeV at the value 
= 3.3 GeV^. It is possible to reproduce the experimental data with fine tuning the continuum 
threshold parameters. However, if we include the contributions of the condensates {qq){q'ggCrGq'), 
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Figure 3: The masses of the scalar mesons as pure tetraquark states with variations of the Borel 
parameter M^, where the (I) and (II) denote the contributions of the condensates {qq){q'gsCrGq') 
of dimension 8 are excluded and included, respectively, q, q' = u, d, s. 


the predicted masses mg are amplified greatly. The mg decrease monotonously and quickly with 
increase of the Borel parameters below some special values, for example, < 1.2 GeV^ 
for the /o(980) and ao(980), then decrease slowly and reach the values ms > 1.4GeV at the 
value = 3.3 GeV^. It is impossible to reproduce the experimental data by fine tuning the 
continuum threshold parameters. In Fig.4, we plot the contributions of different terms in the 
operator product expansion with variations of the Borel parameters for the scalar nonet mesons 
as the pure tetraquark states. From the figure, we can see that the convergent behavior of the 
operator product expansion is very bad, for example, the condensates {qq){q'gsoGq') of dimension 
8 with q, q' = u,d,s have too large negative values at the region > 1.2 GeV^. From Figs.3-4, 
we can draw the conclusion tentatively that the condeiiseLtes{qq){q'gsCrGq') of dimension 8 play an 
important role. The conclusion is compatible with the observation of Ref. [18], that there exists no 
evidence of the couplings of the tetraquark states to the pure light scalar nonet mesons [Si- 

Now we set the mixing angles 0g to be 90° in the QCD spectral densities pg(s) in Eq.(22), 
and take the scalar nonet mesons to be pure two-quark states. In Fig.5, we plot the masses of 
the scalar mesons as pure two-quark states with variations of the Borel parameters the same 
parameters as that in Fig.3 are taken. From the figure, we can see that the predicted masses 
mg ~ (0.85 — 1.14) GeV at the value = (0.5 — 3.3) GeV^, there also exist some difficulty to 
reproduce the experimental data approximately by fine tuning the continuum threshold parameters. 
In Fig.6, we plot the contributions of different terms in the operator product expansion with 
variations of the Borel parameters for the scalar nonet mesons as the pure two-quark states. 
From the figure, we can see that the convergent behavior of the operator product expansion is 
very good, the main contributions come from the perturbative terms, which are of dimension 6 
according to the normalization factors (qq)^ and (ss)^. 

We turn on the mixing angles 6*g ^ 0°, 90° and take into account all the Feynman diagrams 
which contribute to the condensate {qq){q'ggCfGq') with q,q' = u,d,s, see the Feynman diagrams 
in Figs.1-2. The contributions of the vacuum condensates {qq){q'gsCrGq') of dimension 8 can be 
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Figure 4: The contributions of different terms in the operator product expansion with variations 
of the Borel parameter for the scalar nonet mesons as pure tetraquark states, where the 0, 3, 
4, 5, 6, 7, 8, 9 and 10 denotes the dimensions of the vacuum condensates. 



M^(GeV^) 


Figure 5: The masses of the scalar mesons as pure two-quark states with variations of the Borel 
parameter 
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M"(GeV^) 


Figure 6: The contributions of different terms in the operator product expansion with variations 
of the Borel parameter for the scalar nonet mesons as pure two-quark states, where the 6, 
9 and 10 denotes the dimensions of the vacuum condensates. We have taken into account the 
normalization factors (qq)^ and (ss)^. 


13 













































canceled out completely with the ideal mixing angles 0^, 


/o(980) 

= tan“^ 

|0 

ao(980) 

= tan“^ 

o 

o 

00 

o 

O 

= tan“^ 

o 

o 

= tan“^ 


{qq){sgsC^Gs) + {ss){qgsaGq) 

{qq){qgscrGq) 

{qq){sgs<jGs) + {ss){qgsaGq) 

{ss){sgsa'Gs) 

‘^{qq){qgsCrGq) + {qq){sgsaGs) + {ss){qgsCrGq) 


72.6°, 


84.3°, 


{qq){qgsCrGq) 


(4) « 76.0° , 


82.1°, 


( 37 ) 


which results in much better convergent behavior in the operator product expansion. 

In this article, we choose the mixing angles 9s = 0^, then impose the two criteria (i.e. pole 
dominance and convergence of the operator product expansion) of the QCD sum rules on the 
two-quark-tetraquark mixed states, and search for the optimal values of the Borel parameters 
and continuum threshold parameters Sg. The resulting Borel parameters (or Borel windows), 
continuum threshold parameters and pole contributions of the scalar nonet mesons are shown in 
Table 2 explicitly. 

From Table 2, we can see that the upper bound of the pole contributions can reach (51 — 69)%, 
the pole dominance condition is satisfied marginally. If we intend to obtain QCD sum rules for the 
light tetraquark states with the pole contributions larger than 50%, we should resort to multi-pole 
plus continuum states to approximate the phenomenological spectral densities, include at least the 
ground state plus the first excited state, and postpone the continuum threshold parameters Sg to 
much larger values m- In this article, we exclude the contaminations of the continuum states by 
the truncation s^, see Eq.(34), although the truncation Sg cannot lead to the pole contribution 
larger than (or about) 50% in all the Borel windows. Such a situation is in contrary to the hidden- 
charm and hidden-bottom tetraquark states and hidden-charm pentaquark states, where the two 
heavy quarks Q and Q stabilize the four-quark systems qq'QQ and five-quark systems qq'q"QQ, 
and result in QCD sum rules satisfying the pole dominance condition [251 [261127]. 

In Fig.7, we plot the contributions of different terms in the operator product expansion with 
variations of the Borel parameter for the scalar nonet mesons as the two-quark-tetraquark 
mixed states, where the central values of other parameters are taken. From the figure, we can 
see that the dominant contributions come from the vacuum condensates of dimension 6. The 
perturbative contributions of the two-quark components n|^ (p) of the correlation functions IIs (p) 
are proportional to the vacuum condensate (qq)^ (or (ss)^) of dimension 6 according to the nor¬ 
malization factors (qq) (or (ss)) in the interpolating currents Jg(x). In the Borel windows, the 
contributions of the vacuum condensates of dimension 6 are about (109 — 114)%, (90 — 93)%, 
(107— 111)% and (80 — 85)% for the /o(980), ao(980), ko(800) and /o(500), respectively; the con¬ 
tributions of the vacuum condensates of dimension 10 are about (11 — 16)%, (7 — 10)%, (19 — 29)% 
and (16— 22)% for the /o(980), ao(980), ko(800) and /o(500), respectively, where the total contri¬ 
butions are normalized to be 1. The operator product expansion is well convergent in the Borel 
windows shown in Table 2. 

Now we can see that it is reasonable to extract the masses from the QCD sum rules by choosing 
the Borel parameters and continuum threshold parameters shown in Table 2. In Figs.8-9, we 
plot the masses and pole residues of the scalar nonet mesons as the two-quark-tetraquark mixed 
states with variations of the Borel parameters in the Borel windows by taking into account the 
uncertainties of the input parameters. From the figures, we can see that the platforms are very flat, 
the predictions are reliable. In Table 3, we present the masses and pole residues of the scalar nonet 
mesons as the two-quark-tetraquark mixed states, where all uncertainties of the input parameters 
are taken into account. 

There exists a compromise between the minimal masses and the maximal pole contributions, 
and in the following two paragraphs we will show that the mixing angles 9g are optimal values. 
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Figure 7: The contributions of different terms in the operator product expansion with variations of 
the Borel parameter for the scalar nonet mesons as two-quark-tetraquark mixed states, where 
the 0, 3, 4, 5, 6, 7, 8, 9 and 10 denotes the dimensions of the vacuum condensates. 


15 






























































2.0 

1.8 

1.6 

1.4 

C3i,o 


-1-1-1-1-1-1-1-1-1-1-1-1-1-1-r 


2.0 

—1—r 

-1-1-1-1-1-1-1-1-1-1-1-1-r 



^.11 





1 1 1 


. f„( 980 ) 

- Error bound 


1.6 


a„( 980 ) 

- Error bound 

- 

- 


1.4 



- 

- 


CD 1.0 



- 

- 


E 0.8 



- 

- 


0.6 



- 

- 


0.4 



- 

- 


0.2 



- 

. 


0.0 

. 1 

. 



0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 

M^(GeV^) 


0.80 0.85 0.90 


0.95 1.00 1.05 

M^(GeV^) 


1.10 1.15 1.20 


2.0 

1.8 

1.6 

1.4 

S' 1.2 

0 ) 

(3 1.0 
E 0.8 
0.6 
0.4 
0.2 
0.0 


k„(800) 


- Central value 
• Error bound 



0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 

M^(GeV^) 


2.0 
1.8 
1.6 
1.4 

> 1.2 

<D 

C 3 1.0 

E 0.8 
0.6 
0.4 
0.2 


f (500) 


-Central value 

- Error bound 


0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 

M^(GeV^) 


Figure 8: The masses of the scalar nonet mesons as the two-quark-tetraquark mixed states with 
variations of the Borel parameters. 
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Figure 9: The pole residues of the scalar nonet mesons as the two-quark-tetraquark mixed states 
with variations of the Borel parameters. 



Figure 10: The masses of the scalar mesons as two-quark-tetraquark mixed states with variations 
of the mixing angle Os- 
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Figure 11: The pole contributions (PC) of the scalar mesons as two-quark-tetraquark mixed states 
with variations of the mixing angle 9s ■ 



M^(GeV^) 

so(GeW) 

pole 

/o(980) 

0 . 8 - 1.2 

1.5 ±0.1 

(25 - 52)% 

ao(980) 

0 . 8 - 1.2 

1.8 ± 0.1 

(39 - 69)% 

ko(800) 

0 . 6 - 1.0 

1.0 ± 0.1 

(20 - 51)% 

/o(500) 

0 . 6 - 1.0 

1.0 ± 0.1 

(24 - 59)% 


Table 2: The Borel parameters (or Borel windows), continuum threshold parameters and pole 
contributions of the QCD sum rules for the scalar nonet mesons as the two-quark-tetraquark mixed 
states. 


In Fig. 10, we plot the masses of the scalar mesons as the two-quark-tetraquark mixed states with 
variations of the mixing angles 6s, where the input parameters are chosen as = 1.5 GeV^, 

-^/o( 980 ) IGeV , s°Q( 9 go) = 1.8 GeV , Tf^Q( 9 go) = IGeV , = 1.0 GeV , ^^^^(goo) = 

0.8 GeV^, = 1.0 GeV^, = 0.8 GeV^, we introduce the subscripts /o(980), ao(980), 

ko(800) and /o(500) to denote the different Borel parameters. From the figure, we can see that there 
appear minima in the predicted masses at the values ^/o( 980 )/^/o( 98 o) “ 0.6—1.2, ^*ao( 980 )/^a(,( 98 o) = 
0.9-1.1, 6 '«;o( 800 )/^°o( 800 ) = 0.6-1.1, ^'/o( 50 o)/ 6 '°j,( 5 oo) = 0.5-1.2. The lowest masses of the /o(980) 
and ao(980) can reproduce the experimental values approximately; while the lowest masses of the 
ko(800) and /o(500) are larger than the experimental values. In calculations, we observe that the 
minima of the predicted masses vary with the Borel parameters and threshold parameters Sg, 
the mixing angles 9% are the best values. 

In Fig.II, we plot the pole contributions of the scalar mesons as the two-quark-tetraquark 
mixed states with variations of the mixing angles 9s, where the same parameters as that in Fig. 10 
are taken. From the hgure, we can see that the pole contributions increase with the 9s/9^ slowly, 
and reach the maxima at the values 9s/9g = 1.0 — 1.3, then decrease quickly and reach zero 
approximately. The best values appear at the vicinity of the 9g, not far way from the 9g. 

We can draw the conclusion tentatively that the QGD sum rules favor the ideal two-quark- 
tetraquark mixing angles 9g. 

Now we study the finite width effects on the predicted masses. For example, the currents 
Jfg/ao{ 980 )i^) couple potentially with the scattering states KK, we take into account the contri- 
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ms(GeV) 

As(10-4GeV'’) 

/o(980) 


0.98 ±0.06 

8.7 ± 1.3 

ao(980) 


0.97 ±0.05 

5.0 ± 1.7 

«:o(800) 


0.80 ±0.05 

3.6 ±0.6 

/o(500) 


0.70 ±0.06 

5.8 ± 1.0 


Table 3: The masses and pole residues of the scalar nonet mesons as the two-quark-tetraquark 
mixed states. 


buttons of the intermediate KK-loops to the correlation functions nj^g/ap( 98 o)(p^) 
n/o/ao(980)(P^) = -- 


X2 

^/o/ao(980) 


P'' -^%/ao{980)-^KK{p) 


(38) 


where the Xfg/ao( 980 ) Etnd m/o/oo( 980 ) Eire bare quantities to absorb the divergences in the self¬ 
energies All the renormalized self-energies contribute a finite imaginary part to modify 

the dispersion relation, 


X2 

_ ^/o/ao(980) _ _ 

- '”/o/ao(980) + 

The contributions of the other intermediate meson-loops to the correlation functions n 5 (p^) can 
be studied in the same way. 

We can take into account the finite width effects by the following simple replacements of the 
hadronic spectral densities, 


n/o/ao(980)(P^ 


d (s - ml) 


1 v^r5(s) 

7r(s-TO|)Vsr|(s) 


It is easy to obtain the masses. 


ml 


f dss- __ (_ S X 


VsTsis) 


2_ _^ 

^9 ^ (s—mg)^-t-s r^(s) 


exp 


where 


(40) 


(41) 


r/o(980)(8) 

rao(980)('8) 

rKo(800)('S) 

r/o(500)('S) 


r/o(980) , 

Poo (980) ) 

"^ko(800) 
^ ko(800)- 

^ "^/o(500) 

^/o(500)-;- 


(42) 


and the masses ms at the right side of Eq.(41) come from the QCD sum rules in Eq.(34), here 

we have added the factors —and considering the large widths of the ko( 800) and 

/o(500). The numerical results are shown explicitly in Eig.l2. Erom Eig.l2, we can see that the 
predicted masses myQ(g8o) and ^-^^(ggo) are modified slightly after taking into account the small 
widths r/o(g8o) and P^j^i-ggo), the finite widths can be neglected safely; while the predicted masses 
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Figure 12: The masses of the scalar nonet mesons with variations of the Borel parameters after 
taking into account the finite widths. 


m„o(goo) and mj(j( 5 oo) are modified considerably with the largest mass-shifts (5m^(j(8oo) = —0.09 GeV 
and ^TOy(,( 5 oo) = —0.04 GeV. Now the predicted masses from the QCD sum rules are 

m.„( 800 ) = (0.71 ±0.05) GeV, 

"^/o( 5 O 0 ) = (0.66 ±0.06) GeV, (43) 

which are much better than the values presented in Table 3 compared to the experimental data, 

w^o( 80 o) = (682 ± 29) MeV, 

m/„( 5 oo) = (400 - 550)MeV, (44) 

from the Particle Data Group [1]. 

4 Conclusion 

In this article, we assume that the nonet scalar mesons below 1 GeV are the two-quark-tetraquark 
mixed states and study their masses and pole residues using the QGD sum rules. In calculation, 
we take into account the vacuum condensates up to dimension 10 and the 0 {as) corrections to the 
perturbative terms, and neglect the condensates which are vacuum expectations of the operators 
of the order in the operator product expansion. We choose the ideal mixing angles, 

which can lead to good convergent behavior in the operator product expansion, the resulting 
two-quark components are much larger than 50%. Then we impose the two criteria (i.e. pole 
dominance and convergence of the operator product expansion) of the QCD sum rules, search for 
the optimal values of the Borel parameters and continuum threshold parameters, and obtain the 
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masses and pole residues of the nonet scalar mesons. The predicted masses are compatible with the 
experimental data, while the pole residues can be used to study the hadronic coupling constants 
and form-factors. 
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